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ABSTRACT 

The unique properties of dark matter are revealed during collisions between clusters 
of galaxies, like the bullet cluster (IE 0657-56) and baby bullet (MACSJ0025-12). 
These systems provide evidence for an additional, invisible mass in the separation 
between the distribution of their total mass, measured via gravitational lensing, and 
their ordinary 'baryonic' matter, measured via its X-ray emission. Unfortunately, the 
information available from these systems is limited by their rarity. Constraints on the 
properties of dark matter, such as its interaction cross-section, are therefore restricted 
by uncertainties in the individual systems' impact velocity, impact parameter and 
orientation with respect to the line of sight. 

Here we develop a complementary, statistical measurement in which every piece 
of substructure falling into every massive cluster is treated as a bullet. We define 
'bulleticity' as the mean separation between dark matter and ordinary matter, and we 
measure the signal in hydrodynamical simulations. The phase space of substructure 
orbits also exhibits symmetries that provide an equivalent control test. 

Any detection of bulleticity in real data would indicate a difference in the inter- 
action cross-sections of baryonic and dark matter that may rule out hypotheses of 
non-particulate dark matter that are otherwise able to model individual systems. A 
subsequent measurement of bulleticity could constrain the dark matter cross-section. 
Even with conservative estimates, the existing HST archive should yield an inde- 
pendent constraint tighter than that from the bullet cluster. This technique is then 
trivially extendable to, and benefits enormously from larger, future surveys. 

Key vkfords: dark matter — galaxies: clusters: general — cosmology. 



1 INTRODUCTION 

The standard ACDM cosmological model includes a compo- 
nent of cold dark matter that amounts to 85% of the matter 
content of the Universe. This dark matter affects the Uni- 
verse primarily through gravity, and is necessary to explain 
the distribution and growth of large scale structure over cos- 
mic time. However, dark matter does not interact (or only 
very weakly) in the electroweak sector. 

The fundamentally different properties of dark matter 
and baryonic matter are highlighted most dramatically in 
their temporary separation during collisions between galajcy 
cluster pairs, such as IE 0657-56 ([Clowe, Gonzalez fc Marke"^ 



is limited by observational uncertainties in their collision 
velocity, impact parameter and angle with respect to the 
plane of the sky ( Randall et al.| [2008 ) . Current constraints 
are 3 orders of magnitude weaker than constraints from the 



vitch 2004 Clowe et al. 2006 Bradac et al. 20061 and 



MACSJ0025-12 ( [Bradac et al.|2008| ). These 'bullet clusters' 



have provided astrophysical constraints on the interaction 
cross-section a of hypothesised dark matter particles, and 
may ultimate prove the most useful laboratory in which to 
test for any velocity dependence of the cross section. Unfor- 
tunately, the utility of a small number of individual systems 



shapes of haloes (Feng 20101 and, since collisions between 
two massive progenitors are rare ( Shan, Qin fc Zhao|2010 l, 
the total observable number of such systems may be inade- 
quate to investigate a physically interesting regime of dark 
matter properties. 

In this paper, we present a statistical method that al- 
lows every piece of substructure falling into every cluster to 
contribute to a global measure of dark matter-baryonic sepa- 
ration that we refer to as 'bulleticity'. In this approach, every 
infalling mass is treated as a bullet, whose (interacting) gas 
is expected to collide with the intra-cluster medium (ICM) 
and lag behind the (non-interacting) dark matter ( [Powell, | 
Kay fc Babul[[2009 l. Although offsets between baryonic gas 
(e.g. as seen in X-ray emission) and total mass (e.g. as seen 
by gravitational lensing) may not be individually significant, 
detecting a mean bulleticity across many systems would pro- 
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vide robust evidence for a difference between the baryonic 
matter and dark matter interaction cross-sections. Measure- 
ments of the amplitude of buUeticity could then constrain 
the level of the dark matter-dark matter and dark matter- 
baryonic cross-sections. Crucially, since buUeticity should be 
observable in the ongoing assembly of every massive struc- 
ture throughout the Universe, our statistical technique can 
overcome the previous limitations of small number statistics. 
This paper is organised as follows. In Section !|2] we 
develop a simple physical model to illustrate the concept 
of buUeticity and explore some of its dependencies. In Sec- 
tion 33J we use full hydrodynamical calculations to measure 
the expected buUeticity signal in realistic galaxy clusters. In 
Section SE] we discuss the practicality of measurements from 
real astronomical data and use a Fisher matrix analysis to 
predict constraints on the interaction cross-section of dark 
matter using various data sets. We conclude in Section ^ 



2 DYNAMICAL MODEL 

2.1 Motion of test particles falling into clusters 

We shall illustrate the dominant physical effects that sepa- 
rate components of dark and baryonic mass as they fall into 
a massive cluster. Throughout, we shall refer to the com- 
bined infalling system as a 'bullet'. 

Consider the rest frame of a cluster with fixed 'singular 
isothermal sphere' 3D mass distribution p{r) — po{ro/r)^, 
which has a mass interior to radius r of Af (< r) — 4iTporf)r. 
An infalling, point-like component at position r, which in- 
teracts only gravitationally, has an equation of motion 



A-KGporo 



df2 



(1) 



where t is the proper time for the cluster and r is the 3D 
radius. 

An equivalent component that also interacts via the 
electroweak force experiences additional pressure support. 
We approximate this interaction as a buoyancy force equal 
to the weight of the displaced mass. For a sufficiently small 
bullet that the density p{r) of the cluster is constant across 
it, the equation of motion is 



d r 



AnGporo 



1 - (137a) 



2 Poro 



PbT^ 



(2) 



where a is the dimensionless coupling constant and pt, is the 
mean (total) density of the bullet, assuming that the ratio 
of baryons to dark matter is the same in the cluster and the 
bullet. 

We shall use equation ([l| with a = to model the dy- 
namics of standard cold dark matter and equation ([2| with 
a = 1/137 to model baryonic matter. To study interacting 
dark matter, we can simply add a nonzero a term, which rep- 
resents the mean of the dark matter-dark matter and dark 
matter-baryonic coupling constants, weighted by the ratio 
in which the two forms of matter are found in the cluster. 
Such particles would follow orbits between those of stan- 
dard dark matter and baryons. The (137a)^ prefactor can 
be interpreted as either the fractional volume of cluster mass 
displaced by a solid bullet, or as the fractional cross-section 
seen by an interacting particle, i.e. for a geometrically thin 
bullet (137a)^ ^ a/vrr^ 
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Figure 1. Orbits of test particles falling into a 10 Mq 'singular 
isothermal sphere' cluster. In the top panel, the blue curve (large 
arrows) shows the trajectory of standard, non-interacting a = 
cold dark matter, which is governed in this simple model purely 
by the gravitational attraction from the central mass. The red 
curve (small arrows) shows the trajectory of a = 1/137 baryonic 
matter, which experiences an additional buoyancy force. Arrow- 
heads are spaced uniformly in proper time, and the paths pro- 
gressively fade. The bottom panel shows the apparent buUeticity 
if the substructure is moving in the plane of the sky. 



Note that both dark and baryonic matter also experi- 
ence dynamical friction and tidal gravitational forces, but 
equally - so that neither of these effects should separate the 
two components. To keep our model simple, we therefore ig- 
nore these effects, but note that the dissipation of the bullet 
due to tidal and ram-pressure stripping, in conjunction with 
the finite crossing time of a typical cluster, means that we 
would only ever expect to see bullets on their first one or two 
passes through a cluster. Whilst the dynamical friction of a 
lumpy intra-cluster medium would begin to circularise the 
bullet's orbit during this time, we also neglect this effect. 



2.2 Definition of buUeticity 

We define buUeticity b as the vector from the position of the 
dark matter, projected onto the plane of the sky, to that of 
the baryons. Its magnitude is thus sensitive to the difference 
A(T in the total interaction cross-sections of dark matter and 
baryons. Were all matter to have the same interaction cross- 
section, this would ensure b = 0. Any non-zero detection will 
therefore indicate the presence of nonbaryonic material. 

We define the location of the bullet to be the mean 
position (r, 0) of the dark matter and the baryons projected 
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Figure 2. Varying the power law slope of the density profile 
p ex r'^ between 7 = —1 (light) and 7 = —2.6 (dark) in steps of 
0.2, for a fixed cluster mass and bullet density, while adjusting 
the bullet's initial velocity to also ensure a constant impact pa- 
rameter. Notice that db/d'y is small around 7 = — 2 (thick line). 
To avoid confusion, the orbits are only traced up to the second 
periapsis (point of closest approach), and bulleticity to the second 
apoapsis. 



onto the sky. In general, these need not be the same. We 
need to uniquely define the location of each bullet because 
its bulleticity can depend upon the length of time it has been 
within a cluster, and the path it has previously traversed. 
Henceforth, we shall use the symbol r to represent this 2D 
projected radius. 

In the ID case of a bullet that falls into a cluster po- 
tential along a radius, equations (llj) and (pi both produce 
oscillatory motion (with the caveat of a numerical finesse at 
the origin to remove the point of infinite density). However, 
as the baryons experience extra buoyancy, they gradually lag 
behind the dark matter: farther from the cluster core dur- 
ing infall and closer to the core during egress. The absolute 
bulleticity 6 = |b| steadily increases. 

The 2D orbit of a bullet with non-zero impact parame- 
ter but in the plane of the sky is illustrated in Figure [T] In 
this calculation, we assume a cluster of mass 10^'* M© within 
0.8/1"^ Mpc and adopt a value of pb = 10*^ Mq/(/i"^ kpc)^ 
for the baryons; the resulting bulleticity is inversely propor- 
tional to this value. To show the bulleticity more clearly, 
the trajectory is extended to include several passes through 
the cluster, during each of which the separation of dark and 
baryonic matter visibly increases. At each point along the 
mean trajectory, the green line in the bottom panel shows 
the absolute bulleticity b{r). During the first infall, the sepa- 
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Figure 3. Changing the impact parameter from 100 to 
500 h~^ kpc in steps of 50 h~^ kpc (while keeping the impact ve- 
locity and the cluster profile fixed). The default 250 h~^ kpc orbit 
is shown as a thick line. The orbit clearly has a more dramatic 
infiuence than the cluster's density profile (Figurepl, so building 
a theoretical prediction for real clusters will clearly require an 
integration over initial phase space conditions. 



ration gradually increases, with baryons lagging behind dark 
matter as in the ID case. On each subsequent pass, the two 
components separate at large radii, as they follow different 
trajectories, but return to each other as they near the cluster 
core. 

The bulleticity vector can be conveniently written in 
terms of components that are radial br and azimuthal bt 
with respect to the centre of the cluster, such that 



b = bri-r + bti-t 



(3) 



where e^ and et are unit vectors. This is illustrated in figure 
|4] While the radial and tangential bulleticity components 
are smaller than the absolute bulleticity b = |b|, the cru- 
cial point is that they are signed. When averaging over all 
possible viewing angles or the phase space of possible bul- 
let orbits, the mean tangential separation averages to zero 
(it can be easily demonstrated in this case by switching the 
handedness of the orbit to produce a mirror image of Fig- 
ure [IJ. Such symmetries ensure that, for a large sample of 
bullets, (bt) = 0. The remaining, non-zero component of ra- 
dial bulleticity br is shown in Figure IT] as a grey line. 

We thus propose two complementary bulleticity esti- 
mators. The absolute bulleticity b is positive definite, so is 
likely to provide the first detection of a difference in the be- 
haviour of dark and baryonic matter at relatively large S/N. 
The radial bulleticity br has a smaller signal, but subsequent 
measurements from a large cluster sample will benefit from 
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its corresponding statistical check for systematics, bt. This 
latter combination should eventually provide the method's 
cleanest constraints on the interaction cross-section of dark 
matter. 



2.3 Robustness to astrophysical variation 

The distribution of mass in the cluster affects the orbits of a 
bullet. As shown in Figure [2] dark matter precesses around 
a cuspy mass (converging to Keplerian orbits in the limit 
of a point-like central mass, with the mass at one focus of 
the ellipse) , but has constant elliptical orbits around a cored 
cluster, with the centre of the mass at the centre of the el- 
lipse. Importantly however, for a fixed impact parameter, 
the buUeticity of substructure on its first pass through a 
cluster changes by less than 5% for a wide range of clus- 
ter profiles. This lack of influence means that, to first or- 
der, uncertainty about the unknown mass distribution can 
be ignored, and one can freely average results from a large 
ensemble of clusters. It also means that buUeticity measure- 
ments are unlikely to constrain the density profile of galaxy 
clusters. However, other methods are expected to constrain 
this independently ( Massey, Kitching fc Richard|2010 l. 

Conversely, varying a bullet's impact parameter does af- 
fect its subsequent buUeticity. As shown in Figure [3J baryons 
with a low impact parameter pass through more of the intra- 
cluster medium, experience greater buoyancy, and exhibit 
larger buUeticity after periapsis. Bullets with a low impact 
velocity accumulate a greater impulse from buoyancy, so also 
have a higher buUeticity - especially in the tangential direc- 
tion, i.e. bt 3> br- Deriving quantitative predictions about 
the level of buUeticity expected in real clusters will therefore 
require integrations over the phase space of initial conditions 
and subsequent orbits in a typical cosmology, as well as de- 
tailed modelling of the growth and early infall of substruc- 
ture. Such analysis is beyond our simple physical model. To 
do this properly, we shall now switch to full hydrodynami- 
cal calculations, which also automatically include the more 
subtle physical effects that we have disregarded so far. 



3 FULL HYDRODYNAMICAL CALCULATION 

3.1 Properties of the simulations 

We shall now measure the buUeticity signal in realistic, high- 
resolution hydrodynamical simulations of massive clusters 
embedded in a standard cosmological model. The model is 
fiat, with parameters ilm = 1 — ^a ~ 0.3, Qh = 0.04286, 
h = 0.7 and ag = 0.9, where the Hubble constant is defined 
as 100 /ikms~^ Mpc~^ and ag is the mass variance within 
spheres of radius Sh~^ Mpc. All distances are expressed in 
comoving coordinates. 

We performed our simulations with the Adaptive Re- 
finement Tree (ART) A''-body-l-gasdynamics code (|Kravtsov| 



1999 Kravtsov et al. 20021, which uses adaptive refine- 



ment in space and time, and (non-adaptive) refinement in 
mass ( [Klypin et al.|p"001| ) to resolve high dynamic ranges. 
The spatial resolution in the cores of halos is ~ 6/i^^kpc 
and the particle mass is 3-9x 10* h~^ Mq. The formation of 
galaxy clusters is followed from cosmological initial condi- 
tions through later properties of the intracluster medium in- 




Figure 4. Full hydrodynamical simulations of massive clusters 
at redshift z = 0.6. Blue shows total projected mass (dominated 
by dark matter) and red shows X-ray emission from baryonic gas. 
The preferential trailing of gas due to pressure from the ICM, and 
its consequent separation from the non-interacting dark matter, 
is apparent in much of the infalling substructure. 
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Table 1. Mean properties of the simulated clusters. 





'^cluster 


-'Vbullct 


Mass M500 
[1014 Mq] 


Size rsoo 
[h-l Mpc] 


Redshift z = 


= 0.6 








All 


48 


1142 


1.2 ±0.2 


0.75 ±0.05 


Relaxed 


16 


303 


1.6 ±0.5 


0.79 ± 0.09 


Unrelaxed 


32 


839 


1.0 ±0.2 


0.68 ± 0.04 


Redshift z = 


= 








All 


48 


1079 


2.9 ±0.6 


0.73 ±0.06 


Relaxed 


27 


453 


2.0 ±0.6 


0.66 ±0.07 


Unrelaxed 


21 


626 


3.6 ±0.9 


0.79 ±0.08 



eluding gas cooling and star formation (CSF). Astrophysical 
processes in our simulations include metal enrichment and 
thermal feedback due to Type la and Type II supernovae, 
self-consistent advection of metals, metallicity-dependent ra- 
diative cooling and UV heating due to a cosmological ionis- 
ing background. Potentially relevant physical processes ex- 
cluded from the simulations are active galactic nuclei bub- 
bles, magnetic fields, and cosmic rays, although these are 
most important in the innermost cluster regions, which we 
shall exclude anyway. More details about our simulations 



are available in Nagai et al. ( 2007a b) 



We realise simulated 2D observations of 16 independent 
clusters at two redshifts z — 0.6 and z — hy projecting 
snapshots of the dark matter density, gas density and gas 
temperature along three orthogonal axes. This produces 48 
cluster realisations at each redshift, some of which are illus- 
trated in Figure [4] Like the simulations by |Powell, Kay fc| 
Babul (20091, ours show complex interactions of substruc- 



ture. Mock Chandra X-ray imaging was used to visually clas- 



sify each projection as either relaxed or unrelaxed (Nagaietl 
aL][2b07bj . The cluster masses M500 span a range of about 
an order of magnitude centered on lO^'^Af©, and Table n] 
lists the mean mass M500, size rsoo and number of bullets 
A'buiiet in several sub-samples, where the extent of the clus- 
ter defined to be the sphere in which the mean density is 
500 times the critical density of the Universe at that epoch. 



3.2 Analysis of the simulations 

We apply unsharp masking to maps of both the lensing mass 
and the baryonic gas density squared times the square root 
of the temperature (equivalent to the X-ray emission only 
if it is due to thermal bremsstrahlung, but the rest of our 
analysis is completely robust to this simplification). We then 
find the best-fit Gaussian to every local maximum via the 
iteration in size and position adopted by [Rhodes, Refregier| 
|fc Groth] ( |2000[ ). This algorithm is robust but may not be 
optimal in the regime of real observational noise and beam 
smearing. Even in our noiseless simulations, measurements 
near the cluster core are hindered by the steep background 
gradient from central emission, particularly in the X-ray sig- 
nal, which pulls the best-fitting peak inward. More sophis- 
ticated algorithms might prevent this, such as simultane- 
ous fitting of all the sources (c.f. Marshall||2006 l. However, 
we circumvent the issue by ignoring the (very few) peaks 
within 200/i~'^kpc of the cluster centre. We cross-match 
X-ray peaks to their nearest lensing peak projected within 
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Figure 5. Measurements of absolute bulleticity from full hydro- 
dynamical simulations. Error bars show lu errors. Dashed lines 
are best-fitting cubic polynomials to guide the eye. As expected 
for orbits in which only the first crossing has been completed, bul- 
leticity increases as the infalling group falls towards the cluster 
centre and the gas is preferentially retarded. 



6 

a, J 
X 


£^ 

>, -3 
-6 

6 

3 



t -3 

-6 



Q, 



^4 



ff^^ 



w 



0.5 



1.0 1.5 

Projected radius r/r^go 



2.0 



Figure 6. Measurements of radial and tangential bulleticity from 
full hydrodynamical simulations. These quantities are signed so, 
once averaged over many clusters, the values are smaller than 
those in Figure [s] Filled circles show data from redshift 2 = 0.6 
and open circles show data from 2 = 0. The only significant de- 
tection is of radial bulleticity fer(r < rsoo) = 1-7 ± 0.7 h~^ kpc 
at 2 = 0.6. The mean tangential bulleticity is consistent with 
zero, demonstrating its utility as a control test in this harder but 
potentially more discriminating measurement. 



0.5 h,"'^ Mpc, then apply the match the other way around 
and keep only uniquely defined pairs. 

Figures [5] and [6] show the observed mean separation be- 
tween the centres of matched X-ray and lensing peaks, pro- 
jected within 0.5 h~^ Mpc on the sky and averaged over all 
subpeaks. The radial and tangential components br and bt 
can be positive or negative, so they partially cancel upon 
averaging, while b = -^/fop + b^ is positive definite for all bul- 
lets. As in 32J we define the centre of the bullet as the mean 
of the X-ray and lensing positions, and we compute the com- 
ponents of the bulleticity vector with respect to the direction 
towards the global minimum of the projected gravitational 
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potential. As expected from the simple model in 32] the ab- 
solute buUeticity increases towards the centre of the cluster, 
reaching 26/i~^kpc~ 4" at r — O.Srsoo and z = 0.6. The 
normalisation is independent of the cluster mass (at both 
redshifts), within a la uncertainty of 4/i~^ kpc per order of 
magnitude in mass. The measurements are well fit by cubic 
polynomials 



K^) = — 2 ( ^0 + bir + b2r^ + bzr^ 



(4) 



where the coefHcients hi = {38.2,-51.6,35.0,-9.2} a.t z — 
0.6 and h = {24.7, -43.0, 36.4, -10.3} at z = 0, assuming a 
fiducial model in which the prefactor is unity. If the polyno- 
mial order is increased, the quartic coefficients are an order 
of magnitude lower. 

The radial and tangential buUeticity signals are an or- 
der of magnitude lower and noisier. The only statistically 
significant detection is that {br) — 1.7 ± 0.7 h~^ kpc within 
rsoo at 2 = 0.6. This too appears independent of cluster 
mass, and all other measurements are consistent with zero. 
In particular, that the tangential buUeticity signal is always 
consistent with zero demonstrates its utility as a control test 
for any observational systematics. 

3.3 Dependence upon redshift 

The buUeticity signal decreases between the two cosmic 
epochs we have studied. An explanation for this, in agree- 
ment with our predictions from the simple model of ^ is 
that the masses of clusters are considerably lower at z = 0.6. 
With a correspondingly lower infall velocity, the baryons 
therefore experience a larger impulse from buoyancy. Fur- 
thermore, high redshift clusters are also generally more dis- 
turbed, have more substructures per unit mass, and have 
larger offsets between central dark matter and gas peaks. 
However, the projected angle of b on the sky increases for 
nearby clusters. Were the absolute buUeticity to decline lin- 
early with proper time between the values measured at red- 
shifts 2 — 0.6 and 0, the expectation of &(0.3r5oo) would rise 
in angular size to ~ 10" at z = 0.1 and ~ 16" at 2 = 0.05. 

Clusters at low redshift will make the best targets for 
observation. As well as their increasing buUeticity signal, 
nearby clusters also provide a more optimal geometry for 
gravitational lensing. The common misconception to the 
contrary may have arisen from |Hamana, Takada fc Yoshida] 
1 2004 1, who detected peaks using a matched filter of fixed 1' 



size. Nearby clusters appear larger than this on the sky, so 
their fixed filter produces a diminishing signal below z ~ 0.3. 
A better matched filter reveals that the signal is larger, but 
merely spread thin ( Kubo et al.||2007 2009 1 . Finally, mas- 
sive, low redshift clusters should contain more substructure. 
That our sample includes a similar amount of substructure 
at all redshifts is probably a selection effect because our 
high redshift clusters are the most massive progenitors of 
present-day structures, so are growing rapidly. 

3.4 Dependence upon environment 

None of the buUeticity signals show a statistically significant 
dependence upon the cluster's apparent dynamical state. 
However, at redshift z = 0.6 (0), unrelaxed clusters yielded 
about 1.4 (1.8) times as many bullets per cluster as relaxed 



clusters, despite having slightly lower masses at z = 0.6. For 
this measurement, unrelaxed clusters might therefore pro- 
vide more profitable targets. 

We also tried stacking buUeticity measurements from 
different clusters in terms of the bullets' absolute distance 
from the cluster centre, rather than as a fraction of rsoo. 
Because the cluster sizes vary little within our sample, the 
qualitative result does not change. 

In particularly crowded regions of clusters, it would be 
possible to mismatch pairs of dark matter and baryonic pro- 
jected peaks from different substructures. However, this con- 
fusion effect will dampen the measured buUeticity and, if it 
is a function of cluster radius, it will be most pronounced 
towards the centre, where the separation is largest, and the 
density of bullets is highest. If the aim is purely to detect 
buUeticity b in order to prove the existence of dark matter, 
this effect will therefore only make the measurement more 
difficult rather than producing a spurious signal. 



4 PRACTICALITY OF A REAL 
MEASUREMENT 

So far, we have maximally exploited our computationally 
expensive simulations by not adding noise to our mass or 
X-ray maps. We shall now consider the likelihood of and 
practical issues that will be faced by any real measurement 
of buUeticity. 

Baryonic substructure is frequently seen in deep Chan- 
dra imaging (e.g. [Gastaldello et aL||2009| [Randall et al.| 
20091). Dark matter substructure can be mapped efficiently 



via strong gravitational lensing or fiexion (Goldberg & Ba- 
|con|2005| [Bacon et al. 12006] ) , which probes gradients in the 
mass distribution and is therefore more sensitive to small 
mass peaks along a line of sight than weak shear. [Coe et ah] 
( |2010[ ) resolved ten previously unknown subpeaks at various 
radii within (unrelaxed) cluster A1689 using strong lensing. 



and Leonard et al. (2007 20101 resolved four using fiexion 



Both of these measurements were made independently of the 
distribution of light. 

Combined X-ray and gravitational lensing observations 
have already revealed a separation of baryons from dark 
matter in two real systems undergoing major mergers. In 
the bullet cluster IE 0657-56, [Clowe et ah] (J2006') mea- 
sured separations of 6 = 49.3" 
and 46.1" 



Clowe et al. 

: 152 h~ kpc (main cluster) 
142fc~^kpc (bullet). In MACS J0025.4-1222, 



Bradac et al. ( 2008 1 measured separations of dark matter 
from the central gas peak of & = 49.3" = 228 h'^ kpc (SE 
clump) and 30.1" = 139 /i"^ kpc (NW clump). While these 
measurements have small errors, it is their interpretation 
that remains difficult. The extreme disruption of these sys- 
tems has removed any well-defined global potential mini- 
mum, so it is difficult to place the substructure at a well- 
defined radius r or to split the buUeticity into components 
br and bt- Minor merger events will be more usual targets 
for buUeticity. 

Noise in a real measurement - whether due to a finite 
exposure time for X-ray observations or a finite source den- 
sity of lensed background galaxies - will create scatter in 
the measured peak positions. This will emerge as a constant 
minimum b signal. Indeed, this effect is tentatively seen as 
a value around Qh~^ kpc at large r in Figure J5] which coin- 
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cides with the resolution hmit of the simulations. A resolu- 
tion limit will be especially problematic in clusters at high 
redshift, although statistical techniques more sophisticated 
than an offset between peaks will inevitably help. For exam- 
ple, the cross-correlation between the full lensing and X-ray 
maps could be measured. Most importantly, measurements 
of the (signed) radial and tangential bulleticity from a large 
cluster sample will beat down noise on (br), and ensure that 
(bt) -^ 0. 

To predict the observable bulleticity signal, we assume 
that substructure positions can be resolved to 6/i~^kpc, 
which is achievable with Chandra at z < 0.3. At radii where 
the expected bulleticity signal (J4| is resolved, we assume 
that the observed r.m.s. error (7b{r) falls from the base- 
line of Figure pi as l/x/TVft, where A^;, is the number of bul- 
lets. Once the approved Multi-Cycle Treasury programme 
'Through a lens, darkly' (P.I. Postman) has been performed, 
the HST/ACS imaging archive will include ~ 50 clusters be- 
tween redshifts z — 0.1 and 0.2. If these clusters contain just 
~ 100 substructure peaks, and even excluding measurements 
of those within the central 200 h~^ kpc, they should provide 
a detection of b{r) at a signal to noise ratio of ~ 11 when 
integrated over scales 6h~^ kpc < r < 2r5oo- The Euclid sur- 
vey ( Massey et al.|2004 Refregier et al.|2010 l should yield a 
similar detection significance for br{r) and bt{r), even if only 
one substructure peak is identified in each of its ~ 40, 000 
clusters within the same redshift range. 

We can estimate the tightness with which such mea- 
surements will constrain the dark matter interaction cross- 



0.7cm^/g assuming that the main cluster and sub-cluster 



section using the Fisher information matrix ( Tegmark, Tay 



lor fc Hea vens 1997). We again adopt the best fit models 
of b{r) from equation Q, interpolating linearly with proper 
time between z = 0.6 and 2 = 0, and baseline observational 
noise (Ti,(r) around that shown in Figure [s] We assume that 
the bulleticity in radial bins is uncorrelated, and note that 
the dependence on Act is in the mean of 6(r) (not the covari- 
ance). In this case, the one-parameter Fisher information is 



F 



- EE 



bullets r 



dbjr) 
l{r) \dAa 



E (^^2)2 E 



( bo + bir + b2r^ + bar^ )^ 



(5) 



We have summed equally over r bins where r > 200 h~^ kpc 
and b{r) > 6/i~^kpc, but this could be generalised to in- 
corporate a more sophisticated weight function that raises 
the overall signal to noise. For the observational scenarios 
described above, F{'Krl)'^ ~ 100 for the HST archive, and 
- 30, 000 for Euclid. 

Notably, equation ([5| does not depend on the fiducial 
value of Act. With a more comprehensive suite of hydrody- 
namical simulations to more accurately model the behaviour 
of baryonic substructure, we could therefore directly inter- 
pret constraints on Act as those on ct. Assuming typical bul- 
lets of mass mb ~ 5 x 10^^ and radius rt ~ 10 h~^ kpc, we 
thus predict 68% confidence limits of 



1 X 10^25 cm2/GeV for HST 
mb^/F 1 6 X 10"27cm2/GeV for Euclid. 



(6) 



Compare this to 68% confidence limits from the bullet clus- 
ter of a/m < 1.25cm2/g = 2 x 10"24cm2/GeV or cr/m < 



had similar mass-to-light ratios prior to the merger (Ran- 
dall et al.|2008| . 



Even with these fairly conservative estimates, we ex- 
pect a bulleticity analysis of the HST archive to produce 
constraints on ct similar to or tighter than the bullet cluster. 
Furthermore, such constraints are potentially unlimited by 
the uncertainty in orbital parameters for any single object. 
The fundamental strength of our statistical method is the 
trivial way in which it can then be extended to exploit larger 
surveys like Euclid. Bulleticity thus offers a path towards 
ever more discriminating measurements, even if individual 
extreme merger events turn out to be rare ( |Shan, Qin fc| 
Zhao|2010[ ). 



5 DISCUSSION 

We have defined a measure of dark matter-baryonic matter 
separation 'bulleticity' that takes contributions from every 
detected substructure peak in every massive cluster. Any 
non-zero bulleticity measures the difference between the in- 
teraction cross-sections of dark and baryonic matter. This 
interesting new test can be understood via the intuition of 
a simple model, and we have also measured the expected 
value of the bulleticity signal using full hydrodynamical sim- 
ulations of ACDM clusters. 

A conservative estimate of currently available data and 
analysis techniques suggests that there should be enough 
information in clusters from the HST archive to detect 
bulleticity at a signal to noise greater than 10. With further 
hydrodynamical simulations to interpret the absolute level 
of the signal, this could yield independent constraints on 
the interaction cross-section of dark matter at a similar or 
tighter level than the bullet cluster. The real strength of 
this method is the way in which it can subsequently exploit 
large future surveys, free of biases from individual systems 
and in a trivially expandable way. An ambitious, all-sky 
survey could thus rival constraints from particle physics 
experiments (Feng 20101. For a targetted survey, the ideal 



targets would include massive, low redshift, and possibly 
unrelaxed clusters. Most crucially, bulleticity measurements 
will be obtained in a physical regime unapproachable in 
terrestrial laboratories, and may ultimately provide the 
best test for any velocity dependence of the dark matter 
interaction cross-section. 

As a final tantalising prospect, we note that the posi- 
tions of galaxies provide a third (and more easily measured) 



observable. Indeed, Randall et al. ( 2008 1 derived tighter con- 



straints on the dark matter cross section by comparing the 
post-collision locations of dark matter and galaxies (rather 
than dark matter and gas). To first order, galaxies pass 
straight through each other unimpeded because of the sepa- 
ration between them, and between the stars in each galaxy. 
The vector from dark matter to galaxies is therefore a 'bul- 
leticity' measured around a fiducial model of ct = 0, so it 
provides a more direct measurement of nonzero interaction 
cross-section. There are still potential complications to this 



picture. For example, Russell et al. (20101 found a group of 



galaxies in A2146 leading the X-ray emission as expected, 
but the brightest cluster galaxy lags behind it. Complex 
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baryonic physics can also affect observations: [Cortese et al.| 
( 2007 1 showed that stars lead gas in galaxies that are merg- 
ing into clusters, but that new star formation can also be 
triggered in the gas, with the new stars only gradually falling 
forward into the main galaxy. If such effects can be theoreti- 
cally modelled, the most practical tool for this measurement 
is likely to be strong gravitational lensing. If the positions 
of peaks in a standard LensFit mass reconstruction (e.g. 



Richard et al. 2010 1 were allowed to float instead of being 
tied to the positions of galaxies, their measured offsets would 
be precisely this new buUeticity. 
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